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Abstract 

Generalizing the motion representation we introduced for number-conserving 
rules, we give a systematic way to construct a generalized motion represen- 
tation valid for non- conservative rules using the expression of the current, 
which appears in the discrete version of the continuity equation, completed 
by the discrete analogue of the source term. This new representation is gen- 
eral, but not unique, and can be used to represent, in a more visual way, 
any one- dimensional cellular automaton rule. A few illustrative examples 
are presented. 

1 Introduction 

The cellular automaton (CA) model of highway car traffic proposed by Nagel 
and Schreckenberg pQ uses a simple probabilistic cellular automaton rule to 
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represent the motion of cars on a circular one-lane highway with neither 
entries or exits. Despite its simplicity, this model possesses quite a few of the 
essential characteristics of what is observed in a real highway traffic flow. This 
publication was rapidly followed by a huge number of papers on vehicular 
traffic (for a review see [2]). In order to find all the possible highway CA traffic 
rules, we first published a paper [5] in which we established a necessary and 
sufficient condition for a two-state CA rule to be number- conserving later 
followed by a second one [I] in which was proved a more general theorem 
for g-state CA rulesQ Since the direction of motion described by number- 
conserving CA rules may vary in time, number-conserving cellular automata 
(CAs) are more properly models of closed systems of interacting particles. 
When the number of states and number of inputs are not very small, the 
dynamics of the particles is, however, far from being clearly exhibited by the 
rule table. A simpler and more visual picture of the rule can be given by 
its motion representation. The purpose of this paper is to define and give a 
systematic way to obtain such a representation we introduced for number- 
conserving CA rules [3]. 



2 Number-conserving CA rules 

A one-dimensional CA is a discrete dynamical system, which may be defined 
as follows. Let s:ZxNh > Q be a function satisfying the equation 

s(i,t + 1) = f(s(i - r e ,t),s(i -r e + l,t), ...,s(i + r r ,t)), (1) 

for all i £ Z and all i G N, where Z, denotes the set of all integers, N the 
set of nonnegative integers, and Q a finite set of states, usually equal to 
{0, 1, 2, . . . , q — 1}. s(i,t) represents the state of site i at time t, and the 
mapping / : Q r £+ rr+1 — » Q is the CA evolution rule. The positive integers ri 
and r r are, respectively, the left and right radii of the rule. In what follows, 
/ will be referred to as either an (r^, r r ) -radii CA rule or an n-input CA rule, 
where n is the number + r r + 1 of arguments of /. Following Wolfram [6], 

*A possible generalization of this result would be to find a necessary and sufficient 
condition for CAs rule to be eventually number- conserving. Although it is undecidable 
whether a one-dimensional cellular automaton obeys a given conservation law over its 
limit set, a variety of sufficient conditions to be satisfied by a one-dimensional CA to be 
eventually number-conserving have been recently established [5]. 
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to each rule / we assign a code number N(f) such that 



N(f) = J2 X2 ' • • • ' X^- 1 *^-' 3 **™^**. 

When discussing family of specific CA rules it is usually convenient to 
divide the set of all one-dimensional g-state n-input CA rules into equiva- 
lence classes of rules having similar properties. These equivalence classes are 
defined with respect to the group generated by the operators of reflection 
and conjugation denoted, respectively by R and C, and defined on the set of 
all one- dimensional g-state n-input CA rules / by 

, X 2 ; • • ■ ; X n 

) = /(* 

ni Xn— 1) • • • , X\ ) , 

C f{x u x 2 , ■ ■ ■ , x n ) = f(q - 1 - xi, q - 1 - x 2 , . . . , q - 1 - x n ). 

Definition A one- dimensional q-state n-input CA rule f is number-con- 
serving if, for all cyclic configurations of length L > n, it satisfies 

f{Xi,X2, • • • , X n -\, X n ) + f{X2, X3, . . . , X n , X n+ x) + ■ ■ • 

+ f(x L , Xx . . . , X n _ 2 , = Xi + x 2 H h x L . (2) 



In [1] we proved the following theorem. 

Theorem A one- dimensional q-state n-input CA rule f is number-con- 
serving if and only if, for all (x\, x 2 , ■ ■ ■ , x n ) G Q n , it satisfies 



n-1 



f(x lt x 2 ,...,x n )=x 1 +) j (/(0, 0, . . . , 0, x 2 , x 3 , . . . , x n _ k+ i) 



k=l 



-/(0, 0, . . . , 0, x u x 2 , . . . , x n _ k )) , (3) 



For number-conserving CA rules, the motion representation is defined as 
follows. List all the neighborhoods of an occupied site represented by its 
site value s 6 Q. Then, for each neighborhood, indicate the displacements 
of the s particles by arrow(s) joining the site where the particle(s) is (are) 
initially located to its (their) final position(s). A number above the arrow 
indicates how many particles are moving to the final position. To simplify 
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the representation, we only list neighborhoods for which, at least one particle 
is moving. Here are two examples. The motion representation of the two- 
state three-input rule 184, which represents the simplest car moving rule on 
a one-lane highway, is 

i 

10. 

Since a car located behind, that is, on the left of an occupied site cannot 
move, we do not mention it. In a less compact notation we could mention 
this rule completing the representation above with 

o 

11. 

The motion representation of the three-state three-input rule 6171534259461 
is 

112 1 

r\ rx r\ rx 

10 11 20 21 

In both examples particles move only to the right. It is, therefore, not nec- 
essary to indicate the state of the left neighboring site of the particle(s). 

In [7] Hattori and Takesue showed that a CA rule / is number- conserving 
if, and only if, for all {xi,x 2 , . . . ,x n } G Q n it satisfies 

f(x 1 ,X 2 , ...,X n )~ X_ re+1 = J(x U X 2 , ■ ■ . ,X n _i) - J(X 2 ,X 3 , . . .,X n ), (4) 

where 

n— 1 re 

J(x 1 ,x 2 ,..., = —y~] f(0, 0, . . . , 0, xi, x 2 , . . . , x n _ k ) + ^ x r ( 5 ) 

k=l ^ I ' j=l 

Equation (jl]) is actually a discrete version of the continuity equation dp/dt = 
—dJ/dx, where p is the density of the conserved quantity and J its current. 
The left-hand side of (jl]) is the discrete time derivative of the local particle 
density and its right-hand side represents minus the space derivative of the 
current J. More details regarding the current J and its properties can be 
found in [HI [H], and for an extensive discussion of number- conserving rules 
refer to [TT1 H21 E3] ■ 

Here we shall just show how the current J can be used to construct the 
motion representation of number-conserving rules. 

For example, for Q = {0, 1} and r £ = r r = 1, apart the identity, there 
exist four number-conserving rules whose code numbers are 170 (left shift), 
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240 (right shift, reflected of 170), 184 (mentioned above) and 226 (reflected 
and conjugate of 184). For rule 184, the current Ji§4 is given by 

■WO, 0) = 0, J 184 (0, 1) = 0, J 184 (l, 0) = 1, J 184 (l, 1) = 0, 

which can be simply written as 

Jisa{xi,x 2 ) = a?i(l - x 2 ). 

The current is non-zero only for (xi,x 2 ) = (1,0), in agreement with the 
motion representation of rule 184 given above, namely, 

rv 

10. 

The motion representation of rule 226, reflected of 184, is therefore 

01. 

In the case of rule 170, we have 

J 170 (0, 0) = 0, J 170 (0, 1) = -1, J 170 (l, 0) = 0, J 170 (l, 1) = -1, 
or simply 

^170(^15 £2) = ~X 2 , 

and the corresponding motion representation is 

• 1. 

The motion representation of rule 240, reflected of 170, is therefore 

rv 
1«. 

3 Non-conservative CA rules 

For non-conservative rules, (jl]) does not hold, but keeping the definition of 
current J given in ([5]) , we have to add a source term to the right-hand side of 
(HI) to take into account creations and annihilations of particles. If we define 
this source term by 

S(xi, X 2 , ■ ■ ■ , X n ) —f[Xl, X2, • • • , X n ) — X- rt+ i 

, 3?2, • • • , X n — 1 ) - J(x 2 ,x 3 , . . .,x n )), 
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the generalized continuity equation becomes 



f{Xi, X 2 , • • • , X n ) — X- ri+ \ — J [Xi, X2, ■ ■ ■ , X n -i) — J [X 2 , X 3 , . . . , X n ) 

+ Y,(x 1 ,x 2 ,...,x n ). (6) 

Depending upon its sign, positive or negative, the source term S(xi, x 2 , . . . , x n ) 
represents either creation or annihilation of particles. 

In order to illustrate how current and source terms can be used to con- 
struct the motion representation, let us consider again the case of elementary 
CA rules with re — r r — 1. The generalized continuity equation reads 

f(x 1 ,X2,x 3 ) - x 2 = J(xi,x 2 ) - J(x 2 ,x 3 ) + T,(xi,x 2 ,x 3 ), (7) 

where 

J(x!,x 2 ) = -f(0,xi,x 2 ) - /(0,0,xi) +xi. (8) 
For example, for monotone nonincreasing rule 56 we find 

J 56 (0,0)=0, J 56 (0,1) = 0, J 56 (1,0) = 1, J 56 (1,1) = 0, 
S 56 (0,0,0) = 0, S 56 (0,0,1) = 0, S 56 (0,10) = 0, S 56 (0, 1, 1) = 0, 
S 56 (1,0,0) = 0, S 56 (1,0,1) = 0, S 56 (1,1,0) = 0, S 56 (l,l,l) = -1 

The current is the same as for rule 184, but, since E 56 (l, 1, 1) = —1, the 
central particle is annihilated when (xi,x 2 ,x 3 ) = (1, 1, 1). We translate this 
fact in the motion representation by 

r\ —1 

10, 111,. 

In all other cases the particle remains at the same location. The full version 
of the motion representation of rule 56 is therefore 

rv -1 

• 10, 111, Oil. (9) 

To describe the general procedure to construct the motion representation 
from the expressions of J and E we consider the slightly more complicated 
case of rule 174. For this rule, the current is given by 

J m (0, 0) = 0, J m (0, 1) = -1, J m (l, 0) = -1, J m (l, 1) = -1, 
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and the source term by 



£i74(0,0, 0) = 0, £ 174 (0,0,1) = 0, £ 174 (0,1,0) = 0, £ 174 (0,1,1) = 0, 
E 174 (1,0,0) = 1, E 174 (1,0,1) = 1, £ 174 (1,1,0) = -1, £ 174 (1,1,1) = 0. 

Using the current definition, we create the motion representation by orienting 
the arrow either to the right, when the current is equal to 1, or to the left, when 
it is —1. If the current is equal to 0, we omit that entry. For rule 174, this 
yields 

■f\ •Ti 

01, 10, 11. 

The middle term looks rather strange because a particle moves from an empty 
site to an occupied site, but let us not worry about this now since the rep- 
resentation above is temporary. Taking into account the source term this 
anomaly will disappear. Note that the representation above can be written 
in a more compact form as 

• 1, 10. 

Now consider the creation-annihilation part. If £i 74 = 1 (resp. £i 74 = —1), a 
particle is created (resp. annihilated) at the central site and this fact will be 
described writing +1 (resp. —1) above this site. Entries for which S 174 = 
are omitted. Hence, 

E m (l,0,0) = 1 yields lOO, 
£174(1,0,1) = 1 yields 101, 
E 174 (1,1,0) = -1 yields IK), 

Again, this can be written in a more compact form as 

+1 -1 
10«, 110. 

Combining entries obtained from the expressions of Ji 74 and Ei 74 we finally 
obtain 

+1 —l 
• 1, 10, 10«, 110. 

The third term implies the creation of a particle at any empty site located on 
the right of an occupied site, but, according to the second term, this particle 
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will immediately move to its left neighboring site, the motion representation 
of rule 174 can be written in a more compact way 

^ +1 -l 
• 1, 10, 110. 

This representation has to be interpreted carefully. The second term seems to 
imply that a particle is created at an already occupied site if it has an empty 
left neighboring site. This would result in a site occupied by two particles. 
It is, however, not the case since all particles move to the left, and the final 
result of the first and second terms is the creation of a particle at the first 
empty site located at the left end of a sequence of occupied sites. That is, 

^ +i +i 
• 1 and 1 is equivalent to 01. 

We have thus obtained, for rule 174, the motion representation 

+i -i 
01, 110. 

Remark As illustrated by the considerations above, and mentioned in [10] , 
for non-conservative rules, motion representations are not unique. 

4 Conclusion 

Using the notion of current, whose discrete divergence is related to the varia- 
tion of the particle density in the discrete version of the continuity equation, 
completed by the definition of a source term, representing possible creations 
or annihilations of particles, we have shown how to derive from their expres- 
sions a generalized form of the motion representation, we introduced a few 
years ago to describe the temporal evolution of number- conserving cellular 
automaton rules. This new representation, which is not unique, shows in 
a more explicit way how, for each particular neighborhood, particles move 
or are either created or annihilated. A few illustrative examples have been 
given. 
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